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Abstract 

We numerically investigate dynamical and chemical properties of star clusters (open and globular clus- 
ters, and "super star clusters", SSC) formed in interacting/merging galaxies. The investigation is two-fold: 
(1) large-scale (lOOpc-lOOkpc) SPH simulations on density and temperature evolution of gas in interact- 
ing/merging galaxies and (2) small-scale ones on the effects of high gas pressure of ISM on the evolution of 
molecular clouds. We find that the pressure of ISM in merging galaxies can become higher than the internal 
pressure of GMCs (~ 10^ ke K cm^'^), in particular, in the tidal tails or the central regions of mergers. We 
also find that GMCs can collapse to form super star clusters within an order of 10^ yr owing to the strong 
compression by the high presser ISM in mergers. Keywords: Galaxy: kinematics and dynamics 



1 Introduction 

Star clusters (hereafter SC) such as open, globular, and 
super-star clusters are observed to be relatively "minor" 
components (in mass) in luminous galaxies (e.g., the mass 
ratio of the Galactic GCs to the Galactic luminous disk 
mass is an order of 0.1 %). Therefore, it is generally consid- 
ered that dynamical evolution of star clusters (e.g., merg- 
ing or tidal destruction of star clusters) does not infiuence 
significantly the evolution of luminous galaxies. However, 
recent numerical simulations on dynamical evolution of SCs 
have strongly suggested that this SC dynamics is very im- 
portant for the formation of galactic nuclei in dwarf irreg- 
ulars/ellipticals and that of ultra-compact dwarf galaxies 
(Drinkwater et al. 2003), because (1) time scale of dynam- 
ical friction of SCs is relatively short (~ 10^ yr) and (2) 
cross section of SC merging is also large owing to the larger 
Sn of dwarf galaxies (Bekki et al. 2003). Furthermore 
recent HST observations have revealed a significant num- 
ber of young super star clusters (hereafter SSC) in dwarf 
irregular galaxies and suggested that physical properties 
of these SSCs can be correlated with those of their host 
galaxies (e.g., Billett et al. 2002). These recent numerical 
and observational results lead the author to investigate (1) 
how SSCs can be formed in galaxies, (2) how these SSCs 
dynamically evolve, and (3) whether or not there can be 
strong physical relationship between structural and chemi- 
cal properties of SSCs and those of their host galaxies. In 
this paper, we discuss these problems based on numerical 
simulations of SSC formation in galaxies. We here focus 
on SSC formation in major mergers, because recent HST 
observations discovered many SSCs in interacting/merging 
galaxies (e.g., Whitmore et al. 1993). 



2 A new two-fold numerical model 

Since it is very hard to resolve SSC formation sites (less 
than pc) and investigate the formation processes in detail in 
large-scale (10-lOOkpc) galaxy simulations, we adopt 
the following novel two-fold method for the numerical in- 
vestigation of SSC formation. Firstly, we investigate time 
evolution of pressure, temperature, and density of warm 



interstellar medium (ISM) with the initial temperature of 
^ 10* K and the density of 1 atom cm~^, based on galaxy- 
scale (larger-scale) SPH simulations on galaxy mergers with 
total particle number (A*') of ~ 10^. Secondly, we investi- 
gate dynamical evolution and star formation histories of 
GMCs embedded by hot (warm), high pressure ISM, based 
on GMC-scale (< lOOpc) SPH simulations (AT ~ 4 x 10*). 
Primary focus here is to investigate how the high pressure 
ISM (expected for galaxy mergers) trigger collapse of GMCs 
and consequently form SSCs within the GMCs. In the sec- 
ond simulations, the data of the ISM surrounding GMCs 
come from the first simulations: The present simulations 
are fully self-consistent and mimic (or corresponds to) the 
simulations with N ^ 10^. 



2.1 GalEixy-scale 

Our numerical methods for modeling chemodynamical evo- 
lution of galaxy mergers have already been described by 
Bekki & Couch (2001) and the details of the adopted TREESPH 
are given by Bekki (1995). We accordingly give only a 
brief review of them here. We construct models of galaxy 
mergers between gas-rich disks by using the model of Fall- 
Efstathiou (1980). The total mass and the size of a progen- 
itor exponential disk are Md and R^, respectively. From 
now on, all the masses and lengths are measured in units of 
Afd and Rd, respectively, unless otherwise specified. Veloc- 
ity and time are measured in units oiv = (GM^/ Rd)^^^ and 
tdyn ~ {R\/ GMdY^'^ , respectively, where G is the gravita- 
tional constant and assumed to be 1.0 in the present study. 
If we adopt Md = 6.0 x 10^° M© and _Rd = 17.5 kpc as 
fiducial values, then v — 1.21 x 10^ km s~^ and tdyn = 
1.41 X 10* yr, respectively. The dark-to-disk halo mass ra- 
tio and the star-to-gas mass ratio are set equal to 4.0 and 
9.0, respectively. Bulge component is not included in the 
present study. An isothermal equation of state is used for 
the gas with a temperature of 7.3 x 10^ K (corresponding 
to a sound speed of 10 km s"'^). We present the results 
mostly for merger models with nearly prograde-prograde 
orbital configurations, parabolic encounters, and a pericen- 
tric distance of 0.5 -Rd- 
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2.2 GMC-scale 



We numerically investigate the hydrodynamical effects of 
the high pressure ISM on a self-gravitating molecular gas 
cloud in major mergers. The cloud is represented by 20,000 
SPH particles and the initial cloud mass (Md) and size 
(rd) are set to be 10® Mq and 97 pc, respectively, which are 
consistent with the mass-size relation observed by Larson 
(1981). The cloud is assumed to have an isothermal radial 
density profile with p(r) oc l/(r -|- o)^, where a is the core 
radius of the cloud and set to be 0.2rci. An isothermal equa- 
tion of state with a sound speed of Cs is used for the gas, and 
Cs is set to be 4kms^^ (consistent with the prediction from 
the virial theorem) for models with Md = 10® Mq . Each 
SPH particle is also subject to the hydrodynamical force 
of the ISM whose strength depends on the parameters of 
the ISM. The ISM is represented by - 20000 SPH particles 
and and an isothermal equation of state with pressure p (or 
sound velocity of Ch), and density, p. The ISM particles are 
uniformly distributed within a box with size of GRd- A gas 
particle in a cloud is converted into a coUisionless stellar 
particle if (l)the local dynamical time scale [correspond- 
ing to (47rGpi)~'''^, where G and pi are the gravitational 
constant and the density of the gas particle, respectively] 
is shorter than the sound crossing time (corresponding to 
hi/cs, where hi is the smoothing length of the gas), and 
(2) the gaseous flow is converging. This method thus mim- 
ics star formation due to Jeans instability in gas clouds. In 
the model without ISM effects, star formation does not oc- 
cur at all in this star formation method. We mainly show 
the result of the "fiducial model" withp = 1 atom cm^'* and 
P = IQp ke K cm~^, because this model shows the typical 
behavior of ISM-induced star formation in a gas cloud. 



3 Result I. SSC formation sites. 

Figure 1 clearly demonstrates how the gEiseous pressure of 
the interstellar medium becomes dramatically higher (IC*- 
10^ fes) in some regions of the major merger model at T 
= 0.28 Gyr (starburst phase with the star formation rate 
of ~ 75 Mq yr^ ) in comparison with the initial disks. 
Furthermore, gEiseous pressure for some fraction of the g£is 
particles (in particular, in the central regions) exceeds the 
threshold value of 2.0 x 10^ fee for SSC formation, essentially 
because rapid transfer of gas to the central region of the 
merger and efficient gaseous shock dissipation increase the 
gaseous density and pressure greatly. These results clearly 
explain the reason why SSCs can be formed in starbursting 
galaxies and imply that the formation efficiency of SSCs are 
much higher in merging galaxies with massive starbursts 
than in "normal" disk galaxies. Figure 1 also shows that 
SSCs are formed not only in the central region, where gas 
fueling is more efficient, but also in the outer tidal tails and 
"bridges" among the two cores, where the tidal force of the 
merger forms thin and high-density gaseous layers owing to 
gaseous shock dissipation (More details are given in Bekki 
& Couch 2001). 
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Figure 1: Mass distribution projected onto the x-y 
plane in the merger model with at T = 0.28 Gyr for all 
components (upper) and for formation sites (gas or gas 
forming SSCs) of SSCs with P > 10^ kn (lower). In 
total, 2667 SSC formation sites with Pg > 10^ fee (A;b 
is Boltzmann's constant) have been identified prior to 
this starburst epoch. The scale is given in our units 
(17.5 kpc), and each of the frames measures 70 kpc on 
a side. 



4 Result II. Implosive SSC formation. 



Figures 2 and 3 describe how a burst of star formation 
within a gas cloud can be triggered by the external gas pres- 
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Figure 2: Time evolution of a GMC embedded by high 
pressure ISM for gaseous components (cyan) and for 
new stars formed from the gas (magenta). One frame 
measures 200 pc on a side. Note that a very compact 
SSC is formed owing to strong external compression of 
the GMC by the hot, high pressure ISM. 
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Figure 3: Time evolution of star formation rate of the 
model in Figure 2. Note that the time scale of starburst 
is less than a few Myr, which is shorter than the lifetime 
of very massive stars. This implies that supernovae 
feedback effects are less likely to influence super star 
cluster formation. 



sure of the ISM in the fiducial model. In this model with 
moderately strong ISM pressure, the high pressure of the 
ISM can continue to strongly compress the cloud without 
losing a significant amount of gas from the cloud. As the 
strong compression proceeds, the internal density /pressure 
of the cloud can rise significantly. However, the self grav- 
itational force of the cloud also becomes stronger because 
the cloud becomes progressively more compact during the 
ISM's compression. Therefore, the internal gaseous pres- 
sure of the cloud alone becomes unable to support itself 
against the combined effect of the external pressure from 
the ISM and the stronger self-gravitational force. As a re- 
sult of this, the cloud's collapse initially induced by the high 
external pressure from the ISM can continue in a runaway 
manner. 

Due to the rapid, dissipative collapse, the gaseous den- 
sity of the cloud dramatically rises and consequently star 
formation begins in the central regions of the cloud. The 
star formation rate increase significantly to 1.5 Mq yr~^ 
(8 Myr after the start of the cloud's collapse). About 80 
% of the gas is converted into stars within 14 Myr to form 
a stellar system. Because of the "implosive" formation of 
stars from strongly compressed gas, the developed stellar 
system is strongly self-gravitating and compact. This re- 
sult implies that high external pressure from the ISM is 
likely to trigger the formation of bound, compact star clus- 
ters rather than unbound, diffuse field stars. The final mass 
distribution of this model is similar to the King profile, but 
too compact (with the effective radius of ~ 2 pc) to be con- 
sistent with the observed effective radius of typical globular 
clusters (See Figure. 4). 
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Figure 4: Final projected surface density distribution 
of a SSC developed in the model in Figure 2 (in our 
units). Since the system has not yet relaxed dynami- 
cally just after the formation of the SSC, the density 
profile shows a hollow around log (R) = -1.7. 
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5 Conclusions. 

We first demonstrated that SSCs can be formed from GMCs 
"in an implosivo manner" for major mergers: Strong ex- 
ternal pressure of ISM in mergers can trigger the collapse 
of GMCs and consequently cause the rapid star formation 
within the GMCs. The star formation efficiency in the 
GMCs under strong external pressure is as high as 60 — 80 
% (in total cloud mass) so that the developed SSCs can 
be considered to be bound star clusters (e.g., Hills 1980). 
The present study suggests that unbound star clusters that 
are doomed to be destroyed by some external processes 
(e.g., galactic tidal field and interaction with GMC) can 
be formed from GMCs embedded by low (normal) pressure 
ISM whereas SSCs can be formed from GMCs embedded 
by very high pressure ISM. The stars initially within un- 
bound star clusters could soon be identified as "field stars" 
after the destruction of the clusters. The time scale of 
SSC formation is very short (a few times 10® yr, which 
is shorter than the lifetime of very massive stars, which 
are responsible for the type II supernovae) so that the en- 
ergy feedback from supernovae can be much less likely to 
infiuence the SSC formation. Although the "implosive for- 
mation scenario" described here can provide a clue to the 
origin of GCs in general, we have not yet reproduced the ob- 
served Fundamental Plane of GCs (in M31 and the Galaxy) 
and mass-size relation of GCs. Therefore, it is our future 
studies to check whether the implosive scenario can explain 
self-consistently the observed scaling relations of GCs and 
chemical properties of GCs. 
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